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ABSTRACT 

GJ1214b is a warm sub-Neptune transiting in front of a nearby M dwarf star. Re¬ 
cent observations indicate the presence of high and thick clouds or haze whose presence 
requires strong atmospheric mixing. In order to understand the transport and distribu¬ 
tion of such clouds/haze, we study the atmospheric circulation and the vertical mixing 
of GJ1214b with a 3D General Circulation Model for cloud-free hydrogen-dominated at¬ 
mospheres (metallicity of 1, 10 and 100 times the solar value) and for a water-dominated 
atmosphere. We analyze the effect of the atmospheric metallicity on the thermal struc¬ 
ture and zonal winds. We also analyze the zonal mean meridional circulation and show 
that it corresponds to an anti-Hadley circulation in most of the atmosphere with up- 
welling at mid-latitude and downwelling at the equator in average. This circulation must 
be present on a large range of synchronously rotating exoplanets with strong impact on 
cloud formation and distribution. Using simple tracers, we show that vertical winds on 
GJ1214b can be strong enough to loft micrometric particles and that the anti-Hadley 
circulation leads to a minimum of tracers at the equator. We find that the strength of 
the vertical mixing increases with metallicity. We derive ID equivalent eddy diffusion 
coefficients and find simple parametrizations from K^z = 7 x 10^ x m^/s for solar 

metallicity to Kzz = 3 x 10^ x m^/s for the 100xsolar metallicity. These values 

should favor an efficient formation of photochemical haze in the upper atmosphere of 
GJ1214b. 

Subject headings: planets and satellites: atmospheres - planets and satellites: individual 
(GJ1214b) 
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1. Introduction 


GJ1214b is a warm sub-Neptune orbiting a nearby M dwarf and discovered by the MEarth sur¬ 
vey ( Charbonneau et al.|2009 ). It is one of the rare low-mass planets whose atmosphere is character- 
izable in transit spectroscopy by current telescopes. GJ1214b is therefore considered as a privileged 
archetype of this new category of planets also called mini-Neptunes. It has a mass of 6.55±0.98 Mq, 
and a radius of 2.68±0.13 i? 0 , giving it a density of around 1.88 g cm“^ ( Gharbonneau et al.||2009 ). 
This low density necessarily implies the presence of a thick atmosphere. Interior models suggest 
that the planet possesses either a dense iron/rock core surrounded by a thick hydrogen/helium-rich 


atmosphere or a water-rich core surrounded by a thick steam atmosphere (Rogers & Seager 2010 


Nettelmann et al. 2011). Nettelmann et al. ( 2011[ ) found that a water-rich planet depleted in hy¬ 
drogen would require a water-to-rock ratio larger than 6:1. They considered such a large ratio as 
unlikely and favored instead an intermediate case with a hydrogen/helium/water-rich atmosphere. 
GJ1214b has a tight orbit (i.e. semi-major axis of 0.014 AU) and is very likely tidally locked, with 
a permanent dayside and nightside. 

Several observations in transit spectroscopy were performed to probe GJ1214b’s atmosphere 
and to break the degeneracy of possible compositions. First observations by Bean et al. (2010) 


revealed a lack of spectral features between 0.78 and 1 ^m, ruling out the hypothesis of a cloud-free 
solar composition. Other observations ( Bean et al.|2011 Grossfield et al. |2011[[de Mooij et al.|20f^ 


Desert et al.||20lT Berta et al.||2012 Fraine et al.||2013 Nascimbeni et ar]|2015 ) confirmed this flat 
spectrum between 0.6 and 4.5 //m while one group found a significant difference in the transit 
depth measured in J-band (~1.25 /rm) and in Ks-band ( 2.15 /im) ( [Groll et al.||201l1 ). Finally, 
precise measurements with the Wide Field Gamera 3 (WFG3) instrument on the Hubble Space 
Telescope (HST) revealed a very flat spectrum between 1.15 and 1.65 ( Kreidberg et al.|2014b ). 
That implied the presence of high cloud/haze diffusing or absorbing the stellar radiation at low 
pressure. The presence of clouds is required even for a high mean-molecular-mass atmosphere with 
a composition dominated by water, methane, carbon monoxide, nitrogen or carbon dioxide. The 
cloud-top pressure in transit spectroscopy should be less than 10 “^ mbar for a solar-like composition 
and 10“^ mbar for a water-dominated atmosphere (Kreidberg et al. 2014b). 


For the conditions on GJ1214b, these clouds/haze could be either condensate clouds of potas¬ 
sium chloride (KGl) and zinc sulfide (ZnS) or photochemical haze produced by the photolysis of 
methane in the upper atmosphere (Miller-Ricci Kempton et al. 2012; Morley et al. 2013). KGl 


and ZnS clouds are supposed to form between 0.1 and 1 bar (Miller-Ricci Kempton et al. 2012). 


Their presence at low pressure would thus require a strong atmospheric circulation lofting cloud 
particles over several scale heights. The persistence of photochemical haze at low pressure would 
also require a sufficiently strong vertical mixing, counteracting particle sedimentation and bringing 
methane to altitudes where photolysis occurs. Therefore, the formation of high and thick cloud/haze 
on GJ1214b is strongly linked to the atmospheric circulation. The latter has been investigated in 


previous studies using General Girculation Models (GGMs) for cloud-free atmospheres (Zalucha 


et al. 2012; Menou 2012; Kataria et al. 2014). Menou (2012[) used gray opacities and showed that 
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the superrotating zonal winds and the thermal phase curves would vary greatly with the atmo¬ 
spheric metallicity. Kataria et al. (2014) analyzed the atmospheric dynamics and thermal phase 
curves with a more accurate GCM using non-gray opacities for different atmospheric compositions, 
but they did not analyze the vertical winds and the vertical mixing. 


In this paper, we analyze the circulation and the vertical mixing of GJ1214b for a cloud- 
free atmosphere with different compositions using the Generic LMDZ GGM, a very versatile 3D 
model developed to simulate any kind of planetary atmosphere. In the next section, we describe 
the model and the parameters used for GJ1214b. In section 3, we analyze the thermal structure 
and atmospheric circulation (zonal, meridional and vertical winds). In section 4, we analyze the 
vertical mixing of tracers and discuss its implications for cloud/haze formation in section 5. The 
3D modeling of realistic clouds with our GGM will be described in a dedicated paper, analyzing 
their formation, dynamics and impacts on spectra. 


Model 


2.1. The Generic LMDZ GCM 


We simulated the atmosphere of GJ1214b using the Generic LMDZ GCM. This model has 
been specifically developed for exoplanet and paleoclimate studies. It has been used for studying 


the early atmospheres of the Earth, Mars and Titan (Charnay et al. 2013 Wordsworth et al 


2013 

Forget et al. 

2013 

Charnay et al. 

2014) and for studying the habitable zone for terrestrial 

exoplanets ( 

Wordsworth et al.||2011 

Leconte et al.||2013a|b). The model is derived from the LMDZ 

Earth GGM ( 

Hourdin et al.||2006 

), which solves the primitive hydrostatic equations of meteorology 


using a finite difference dynamical core on an Arakawa C grid. This dynamical core has been used 


in GCMs dedicated to the present Earth (Hourdin et al. 2006), Mars (Forget et al. 1999), Venus 


(Lebonnois et al)]|2010), and Titan (Lebonnois et al. 2012). 


In this paper, simulations were performed with a horizontal resolution of 64x48 (corresponding 
to resolutions of 3.75° latitude by 5.625° longitude). We also did a few tests for GJ1214b with a 
128x96 resolution but we did not notice significant differences. For the vertical discretization, the 
model uses pressure coordinates. In this work, we used 45 layers equally spaced in log pressure, 
with the first level at 80 bars and the top level at 3 Pa. 

The radiative scheme is based on the correlated-k model, with k-coefficients calculated from 
high resolution absorption spectra, computed by kspectrum, a line-by-line model using the HITRAN 
2012 database (Rothman et al. 2013). kspectrum is a tool developped by Vincent Eymet and 


available online at http://www.meso-star.com/en_Products.htnil. For water opacities, we used 
the high-temperature database HITEMP 2010 ( Rothman et al.||2010 ). For methane, the HITRAN 
database does not include absorption lines between 9000 and 11000 cm“^ and between 11500 and 
12000 cm“^. We added the missing lines using EXOMOL (Tennyson &: Yurchenkoj2012 Yurchenko 








































&: Tennyson||2014). At a given pressure and temperatures, correlated-k coefficients in the GCM are 


interpolated from a matrix of coefficients stored in a 12 x 9 temperature and log-pressure grid: T = 
100, 200, 1000, 1200, 1500, 2000 K, P = IQ-^ , 10° , 10^ , 10^ Pa. To facilitate comparison 

with the work of Kataria et al. (2013), we computed k-coefficients over the same 11 spectral bands. 
The H 2 -H 2 and H 2 -He CIA from the HITRAN database were included. Rayleigh scattering by H 2 , 
H 2 O and He was included, based on the method described in Hansen &: Travis (1974) and using 
the Toon et al. (1989) scheme to compute the radiative transfer. 


The model also includes a scheme for adiabatic temperature adjustment triggered in the regions 
where the atmosphere is convectively unstable. However, this condition almost never occurs in our 
simulations of GJ1214b. 


2.2. Atmospheric composition 

We ran simulations with four different atmospheric compositions: 

• H 2 -rich atmosphere at lx solar metallicity 

• H 2 -rich atmosphere at 10 x solar metallicity 

• H 2 -rich atmosphere at 100 x solar metallicity 

• pure H 2 O atmosphere 


For the H 2 -rich atmospheric composition, we used the solar nebula atomic abundances from 


Fodders (2003). The 10 x solar metallicity is similar to Saturn’s metallicity, and the 100 x solar is 


quite similar to Uranus and Neptune’s metallicity (Kreidberg et al. 2014a). Population synthesis 


models generally predict high metallicities between 100-400 x solar for low-mass planets but lower 


metallicities can also arise (Fortney et al. 2013) 


For the H 2 -rich atmospheres, we computed the compositions at thermochemical equilibrium, 
assuming only H 2 , He, H 2 O, GH 4 , N 2 , NH 3 , CO and C 02 . We computed the mixing ratios analyt¬ 


ically following the method described in the appendix of Burrows Sz Sharp (1999), using molecular 


Gibbs energies from Sharp Sz Huebner (1990). We computed abundances from the reactions: 


CO + 3 H 2 = H 2 O + CH 4 

( 1 ) 

N 2 + 3 H 2 = 2 NH 3 

( 2 ) 

CO 2 + H 2 = CO + H 2 O 

(3) 


The computation requires solving a cubic equation. 


Figure ?? shows the vertical molecular abundances computed for 1, 10 and 100 times the solar 
composition cases using temperature profiles from the ID model (see figure [^. 
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The opacities for a H 2 -rich atmosphere are dominated by molecular absorption (mainly H 2 O). 
However, at high pressure and for temperatures higher than around 1000 K, the pressure broadening 
of atomic alkalin (Na and K) spectral lines has a strong impact on the opacity in the visible. For a 
solar composition, atomic potassium should be the dominant K-bearing gas for temperatures higher 
than around 980K at 1 bar ( Fodders Sz Fegleyf2006 ). We did not include Na opacities in our model, 
as Na absorbs at wavelengths that are shorter than the bulk of the radiation emitted by a later-type 
M dwarf like GJ1214. Also, at the high metallicities considered here, Na is likely to condense as 
Na 2 S clouds in the deeper atmosphere ( Morley et al.||2013 ), removing it as an atmospheric opacity 
source. For K, we simply added a continuum between 1.18x10^ - 1.43xl0‘^cm“^, assuming pressure 
broadening and a reference cross section for atomic K of 10“^®cm^/atom at 1 bar and 1000 K for 
solar metallicity. This value corresponds to the absorption from Burrows & Volobuyev (2003) at 
the edges of the spectral band we considered. 


2.3. Integration 


Model parameters are given in tables 1 and 2. For the stellar spectrum, we used a blackbody 
at 3026 K. We also tested with spectra of AD Leo and Gliese 581 but these more realistic spectra 
did not produce significantly different results. We assumed a stellar flux of 23600 W/m^ at the top 
of the atmosphere for a null zenith angle. We also assumed an internal luminosity of 0.73 W/m^ 
corresponding to an intrinsic effective temperature (i.e. the blackbody temperature corresponding 
to the internal luminosity) of 60 K, as suggested by Rogers & Seager (2010). 


As initial states for the 3D simulations, we used the ID temperature profiles computed with 
the ID version of the model (see figure and with no wind. The ID model uses the same vertical 
discretization and physics (e.g. radiative transfer, convective adjustment) as the GCM, assuming 
a stellar zenith angle of 60° and a redistribution of energy over the entire planet. 


For the ID calculation, we ran the model for 1000 GJ1214b days (i.e. 1600 days) and increased 
the radiative heating/cooling by a factor proportional to the pressure when the pressure was higher 
than 0.1 bar. This enables faster convergence and is equivalent to running the model for around 
10^ days for the deepest level at 80 bars. This technique resulted in simulations very close to 
equilibrium after these 1600 days of integration, with relative differences between total emitted 
radiation and total absorbed radiation lower than 0.1% and with no observable variation in the 
temperature profile for longer simulations. The 3D simulations were also run for 1600 days but we 
did not use the technique to accelerate the convergence since the atmosphere is not in radiative 
equilibrium in 3D. The differences between total emitted radiation and total absorbed radiation for 
the final 3D states were lower than 1%. 


Simulations were performed with a dynamical timestep of 60 seconds and a physical/radiative 
timestep of 300 seconds. The GGM outputs are either instantaneous values or daily averaged values. 
All the results we present here use daily averaged values. However, the 3D simulations show small 










time variability with, for instance, standard deviations for temperature generally lower than 10 K. 


3. Atmospheric dynamics 
3.1. Thermal structure 


The left panel in figureshows the temperature profiles obtained from the ID model to initial¬ 
ize the 3D runs. Dashed and dotted lines correspond to the saturation vapor pressure curves for KCl 


and ZnS for the solar and lOOxsolar metallicity from Morley et ah (2012). The temperature prohle 


is only adiabatic below the isothermal region (i.e. below 40 bars for the lOOxsolar composition and 
below 100 bars for the other cases) where it is controlled by the internal heat flux. For these ID 
temperature profiles, KCl and ZnS should condense at around 0.4 bar for the solar metallicity and 
at around 60 mbar for the lOOxsolar metallicity. These altitudes are similar to those obtained by 


previous ID models (Miller-Ricci Kempton et al. 2012; Morley et al. 2013) and are fairly strongly 


dependent on the metallicity. Our 3D simulations globally give the same condensation altitudes, 
however there are latitudinal variations because of the equator-pole temperature gradient. For 
instance, for the lOOxsolar metallicity (the case with the strongest variations), KCl condensation 
would start at around 40 mbar at the equator and at around 200 mbar at the poles. 

In the 3D simulations, there are strong day-night temperature variations that become more 
pronounced with increasing atmospheric metallicity (see right panel in figure [^. At 0.1 mbar the 
day-night temperature difference reaches 500 K for the water atmosphere but is lower than 100 
K for the solar metallicity. The atmospheric opacity increases with metallicity while the specific 
heat capacity decreases. Both effects lead to larger heating/cooling rates above 0.1 bar for higher 
atmospheric metallicity (see figureleft). The altitude where stellar energy is deposited increases 
with metallicity (see figure [fright). In our model the pure water and the lOOxsolar atmospheres 
have similar net stellar flux. The lower water vapor absorption in the lOOxsolar atmosphere 
compared to the pure water atmosphere is compensated by absorption by other gases, in particular 


methane. With our model, we obtained net stellar and thermal fluxes similar to those in Kataria 
et al. (2014) for our four different cases. However, the day-night temperature difference for the 


solar metallicity is a factor of two weaker than in their model. This may be primarily due to the 
stronger equatorial jet we obtained for this case (see next paragraph). 


3.2. Zonal winds 


Concerning the zonal winds, all simulations lead to the development of an equatorial superro- 
tating jet as expected for strongly irradiated tidally locked exoplanets ( Showman &: Polvani|[20TT ). 
Figure shows the zonally averaged zonal wind for our four atmospheric compositions. For the 
lx solar metallicity, there are also two high-latitudes jets. These jets are strongly reduced for the 

















10 X solar metallicity and vanish for the 100 x solar metallicity. Their development is controlled by 
latitudinal temperature gradients through the thermal wind equation (Lewis et al. 2010). With 


higher metallicity, stellar energy is absorbed higher in the atmosphere leading to weaker latitudinal 
temperature gradients in the deep atmosphere and so to less significant and shallower jets. In 
contrast, the depth of the equatorial jet tends to increase with the metallicity leading to stronger 
dayside-nightside temperature gradients. This behavior (increase of the equatorial jet and decrease 


of the high-latitude jets with metallicity) was also obtained in the simulations of Kataria et al. 


(2014). 


Compared to the results from Kataria et al. (2014), our GCM produces shallower high-latitude 


jets for the solar and 10 x solar metallicity cases and no jet for the 100 x solar metallicity. In their 
model, the high-latitude jets reach a maximum speed between 10“^ and 10“^ bar, while maximum 
speed occurs at the top of our model (at 0.03 mbar) for the solar and 10 xsolar metallicity case. 
Therefore, high-latitude jets are stronger at low pressures in our model but weaker at high pressures. 
For instance, for the solar metallicity, the high-latitude jet speed is around 1100 m/s and 500 m/s 
at respectively 10“^ and 10“^ bar in our model whereas it is around 900 m/s and 600 m/s at 


the same pressures in Kataria et al. (2014). Concerning the equatorial jet, our model produces 
a stronger equatorial jet with a wind speed of 1500 m/s at 10“^ bar compared to 900 m/s in 


Kataria et al. (2014). However, the equatorial jets are very similar between both models for higher 


metallicity. These differences between both models for H-rich atmospheres may be due to some 
variations in atmospheric opacities. Otherwise, it may be due to the different dynamical cores, 
based on a standard longitude-latitude grid with a small dissipation in the Generic LMDZ and 


on a cubed-sphere grid with no dissipation in Kataria et al. (2014). However, our results globally 
remain similar to theirs. 

Concerning the pure water atmopshere, our model produces an equatorial jet weaker than for 
H-rich atmospheres, with a maximum wind speed of around 1 km/s compared to 2 km/s for the 
others cases. For the water atmosphere, there also are two weak high-latitude jets with a maximum 


speed of around 500 m/s at 70°N/S and at 1 mbar. Kataria et al. (2014) obtained the same 


structure and speeds for the zonal wind in the pure water atmosphere. 


3.3. Meridional circulation and vertical winds 

For a tidally locked, synchronously-rotating exoplanet, the heat redistribution from dayside 
to nightside is mostly done by the superrotating jets. There is also a meridional circulation re¬ 
distributing heat from low to high latitudes. This is illustrated in figure which shows - for the 
100 X solar metallicity - equatorial temperature and vertical wind as a function of pressure and 
longitude, as well as maps of temperature and winds at 1.1 mbar (altitude close to the infrared 
photosphere). The dayside is characterized by upwelling centered near the sub-solar point with 
poleward winds transporting heat to high latitudes. The dayside poleward circulation crosses poles 
and is equatorward on the nightside, allowing a stronger heat redistribution. This circulation is 






















shaped by standing Rossby and Kelvin planetary waves, responsible for the superrotation by pump¬ 
ing eastward momentum from high latitudes to the equator ( Showman &: Polvani|20lT ). Winds and 
temperatures are characterized by a chevron-shape pattern, centered at the equator and pointing 
to the east (Showman & Polvani 2011). Two strong downdrafts are present at the equator close 
to the terminators. They are associated with an adiabatic heating explaining the hot spot west of 
the anti-stellar point in Fig. They are also associated with a shock-like feature with a strong 
updraft at the substellar and anti-stellar points, also observed in 3D simulations of hot Jupiters 
( Showman et al.||2009t Rauscher Sz Menou 2010). 


Figure shows the mass streamfunction for the 100 x solar metallicity case for the dayside, the 
nightside and globally (i.e. dayside-|-nightside). For the dayside, we integrated longitudes between 
-39° and -1-90° E only, to avoid the very strong downdraft located west of the substellar point. 
For the nightside, we integrated the mass streamfunction over all the longitudes not included in 
the dayside integration. The dayside is characterized by a large Hadley circulation for pressures 
below 0.1 bar with two cells extending from the equator to the poles. On the nightside, the mean 
meridional circulation reverses. This reversal is due to equatorward winds on the nightside and to 
the two very strong downdrafts located at the equator and close to terminators. These downdrafts 
are typically twice as strong as the equatorial updrafts (see figure]^. With an integration between 
-90° and -|-90° E for the dayside, so including the downdraft located west of the substellar point, 
we obtained similar circulations but weaker at the equator. When looking at the global mass 
streamfunction, it appears that the anti-Hadley circulation is stronger in the deep atmosphere (i.e. 
between 0.1 and 0.01 bar) than the dayside Hadley circulation, but weaker in the upper atmosphere. 
In reality, the meridional circulation corresponds rather to two large cells crossing the poles and 
going from the dayside to the nightside but the division into a dayside Hadley circulation and a 
nightside anti-Hadley circulation allows an evaluation of which one dominates. These results are 
for the 100 X solar atmospheric composition, and yet the other atmospheric compositions present 
the same circulations for pressure lower than around 0.1 bar. 


Figure shows a few contours of the mass streamfunctions averaged zonally over all longitudes 
and superposed on the zonal mean vertical winds. The mass streamfunctions follow variations of 
the zonal mean vertical winds. For the lx and 10xsolar metallicity at pressures lower than 0.1 bar, 
the zonally averaged circulation is characterized by an upwelling at mid-latitudes (between around 
20°-60°) and a downwelling at the equator. The zonal mean circulation is thus dominated by the 
nightside anti-Hadley circulation. The circulations for the 100 x solar metallicity and the pure water 
atmospheres are more complex, with the anti-Hadley circulation dominating at pressures greater 
than 10 mbar and the Hadley circulation above. As detailed before, the dayside Hadley circulation 
becomes stronger than the nightside anti-Hadley circulation in the upper atmosphere, explaining 
this transition. Therefore the zonal mean anti-Hadley circulation is present in every case but it is 
located deeper in the atmosphere for higher metallicity. 


The top panel in figure shows the zonal mean difference between the outgoing longwave 
radiation (OLR) and the absorbed stellar radiation (ASR). Positive values correspond to regions 










_ 9 _ 


with surplus energy (compared with the radiative equilibrium) transported by the atmospheric cir¬ 
culation from regions with a deficit of energy (negative values). Globally, the meridional circulation 
transports heat from low latitudes to high latitudes. However, the anti-Hadley circulation limits 
this poleward transport at low latitudes. The equatorial downdrafts produce an adiabatic heating 
(see for instance the equatorial temperature at 150°E in figure Q. The difference between OLR 
and ASR exhibits a local maximum at the equator due to this adiabatic heating. For the 100 x solar 
metallicity, the meridional circulation contributes to heating at the equator in average. 

The bottom panel in figure [7] shows the vertically and zonally-integrated heat flux and its 
decomposition into components of the heat flux transported by mean meridional circulation, sta¬ 
tionary waves and transient perturbations. We computed the heat flux from the meridional wind 
V and the dry static energy e = (cpT -|- gz), and we decomposed the heat flux as: 

[eu] = [e][u] -|- [e*v*] + [e'][u'] (4) 

where square brackets are time averages and asterisks are deviations from time averages, overbars 
are zonal averages and primes are deviations from zonal averages. The terms on the right represent, 
respectively, the transport by the mean meridional circulation, the transport by the stationary waves 
and the transport by the transient (time-dependent) perturbations. For the 100xsolar metallicity 
case, the total heat transport is strongly impacted by the sationary waves. At a given pressure level, 
the dry static energy is higher (both T and z are higher) on the dayside with net poleward winds 
than on the nightside with net equatorward winds. Therefore, the stationary waves associated with 
the day-night contrast produce a net poleward heat transport, counterbalancing the equatorward 
heat transport associated with the mean anti-Hadley circulation below 10 mbar. 

Finally, the anti-Hadley circulation we described in this section should be present on a large 
range of warm tidally locked exoplanets and has a fundamental impact on vertical mixing and cloud 
transport. 


4. Vertical mixing 

We performed an analysis of the strength of the vertical mixing in GJ1214b’s atmosphere 
using the Generic LMDZ GGM with simple tracers. The goal was to understand how the general 
circulation, and in particular the zonally averaged circulation described in the previous section, 
controls the transport and distribution of potential cloud or haze particles. From this 3D analysis, 
equivalent ID eddy diffusion coefficients were obtained, which may be of use for ID chemical models. 


4.1. Implementation in the GCM 

We used a similar method to the analysis of vertical mixing in the hot Jupiter HD 209458b by 


Parmentier et al. (2013). We added particles with a fixed diameter and with a constant abundance 
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no below 1 bar, with this region corresponding to the source/sink of particles. We let the circulation 
transport these particles and allowed them to sediment at a terminal velocity described in the next 
section. Contrary to Parmentier et al. (2013), we do not assume that particles evaporate on 


the dayside, which would reduce the downward sedimentation tracer flux. Retaining the dayside 
particles is based on our model dayside temperatures, which are not high enough to significantly 
evaporate KCl or ZnS clouds even for the 100xsolar metallicity case. 

We ran simulations for the different atmospheric compositions with particle radii from 0.1 
to 10 microns and with a particle density corresponding to solid KCl (2000 kg/m^), a possible 


composition for GJ1214b’s clouds (Miller-Ricci Kempton et al. 2012). We considered micrometric 


particles because they correspond to the estimated size for possible high cloud/haze particles that 
would produce a flat transit spectrum in visible and near-infrared (Miller-Ricci Kempton et al.|2012 


Morley et al. 2013). For the initial state, we assumed a uniform particle abundance everywhere in 
the atmosphere and equal to no (the constant abundance below 1 bar). We used no=10“® kg/kg 
in all simulations. This value is arbitrary and the next figures always show the relative abundance 
n/no- 


4.2. Sedimentation 


We computed the sedimentation of particles assuming that they fall at the terminal velocity 


given by (Fuchs 1964; Ackerman & Marley 2001); 


^/ = 


Wr’^gjPp - p) 
9r] 


(5) 


where r is the particle radius, g is the gravitational acceleration of GJ1214b, pp is the particle 
density (we chose pp=2000 kg/m^ corresponding to solid KCl), p is the atmosphere density, p is the 
viscosity of the atmospheric gas and /3 is the Cunnigham slip factor, which describes non-continuum 


effects. An experimental expression of the Cunnigham slip factor is (Fuchs 1964): 


/3 = 1 + Kn (l.256 + 0.4e"^-^/^") 


( 6 ) 


where Kn is the Knudsen number, equal to the ratio of the mean free path to the size of radius 
of the particle: 


The mean free path A is given by: 


K„ 


A 

r 


(7) 


ksT 1 
y/l-KcP P 


( 8 ) 
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with ks the Boltzmann constant, T and P the temperature and pressure of the gas and d the 
molecular diameter. 


For the hydrogen-rich atmospheres, we considered that the dynamical viscosity was equal to 
that of pure H 2 . Following Ackerman &: Marley (2001) and Parmentier et al. (2013), we used the 
formula: 


m2 = 


5 V^rnksT {kBT/e)^-^^ 


(9) 


16 nd'^ 1.22 

where m is the molecular mass, d=2.827x10“^° is the molecular diameter and e=59.7ks K is the 
Lennard-Jones potential for H 2 . The ratio on the right in the formula corresponds to the deviation 
from the hard-sphere model. 

H 2 O is a polar molecule that deviates too strongly from the hard-sphere model to be described 
by formula ( 8 ) (Crifo|19M). For the water atmosphere, we used a parameterization of the dynamical 


viscosity of H 2 O for dilute conditions and given by Sengers & Kamgar-Parsi (1984): 


T] 


^HaO = 




k=0 


( 10 ) 


with 7]*= lO-f^ Pa/s, r*= 647.27 K, ao= 0.0181583, ai= 0.0177624, 02 = 0.0105287 and 03 = 
-0.0036744. This parameterization is valid for temperatures ranging from around 300 K to 1100 K. 
The dilute approximation is valid up to 10 bars for GJ1214b conditions (see Sengers & Kamgar- 


Parsi (1984) for an expression with correction terms to the dilute approximation), and is thus 


compatible with the pressure range we considered for vertical mixing. For the calculation of 
for the water-rich atmosphere, we used a molecular diameter dH 20 depending on temperature and 
given by ( Crifo||1989 ): 


'T' \ — 0.3 

300kJ 

The top panel in figure shows the terminal velocity for particle radii from 0.1 to 10 microns, 
using /?p=2000 kg m“^ (similar to the density of KCl) and for the 100 x solar and the pure water 
atmospheric composition cases (temperature profiles from the ID model). For both compositions, 
two regimes are present. At high pressure, the Knudsen number is smaller than 1. The Cunningham 
slip factor and thus the terminal velocity are almost constant. At low pressure, the Knudsen number 
is much larger than 1. The Cunningham slip factor and thus the terminal velocity are inversely 
proportional to the pressure. For a given particle radius, the terminal velocity and the corresponding 
sedimentation rate are always lower in the H 2 O atmosphere than in the H-rich atmosphere. In the 


dH20 = 4.597 X 10“^°m 

























12 - 


upper atmosphere, the terminal velocities are lower by around a factor 3. This is due to both the 
larger molecular diameter and the larger viscosity of H 2 O compared to H 2 . 

The sedimentation timescale is defined as the ratio of the scale height to the sedimentation 
velocity Ts = H/Vf. The bottom panel in figure shows the sedimentation timescales for the 
100 X solar and the pure water atmosphere. The vertical mixing simulation is valid when the duration 
of the simulation is longer than the sedimentation timescale. Since our simulations were run for 
around 1600 days, the vertical mixing simulations are valid for pressures lower than 0.1 bar for 
a particle radius of 1 micron, and for pressure lower than 1 bar for particle radii of 3-10 micron. 
For evaluating the vertical mixing of the tracers, it is useful to perform simulations with different 
particle radii. Large particles provide the most robust results at higher atmospheric pressures, 
but can be completely removed from the upper atmosphere, and so smaller particles give more 
accurate estimates for vertical mixing at these levels. Figure shows the ratio (in log scale) of 
the sedimentation timescale to the advection timescale {Tadv = 2ttR/U with R GJ1214b’s radius 
and U the mean zonal wind) for the 100xsolar case. This ratio is very similar for the different 
atmospheric compositions. When the sedimentation timescale is longer than the advection timescale 
(i.e. particles fall by less than a scale height during a full rotation around the planet), the particle 
transport is mostly driven by the zonal mean circulation and there are small longitudinal variations 
in the particle distribution. According to figure the particle distribution is mostly controlled by 
the zonal mean circulation and there should be almost no variations in longitude for pressures 
higher than 1 mbar. Moreover, longitudinal variations should be enhanced at high latitude. 


4.3. Results for the 3D mixing 

Figure [To| shows the global mean particle abundance as a function of pressure for the different 
atmospheric compositions for a particle radius of 1 micron and for the 100 x solar metallicity case 
for particle radii from 0.1 to 10 micron. The particle abundance in the upper atmosphere increases 
with atmospheric metallicity. The abundance is reduced by 50% at around 5 mbar for the solar and 
10 Xsolar metallicity cases, and at 0.2-0.4 mbar for the other cases. The strength of the atmospheric 
circulation above 0.1 bar increases with metallicity. For the solar metallicity case, the circulation 
cannot efficiently transport micrometer particles to the upper atmosphere (0.1-0.01 mbar). For 
the 100 Xsolar metallicity, the vertical mixing is weak between 1 bar and 0.1 bar and very strong 
between 0.1 bar and 1 mbar. Particles with a radius of 0.1 micron are well mixed above the particle 
source layer while particles with a radius of 10 microns cannot go higher than 0.5 bar. 

Figure 0 shows the zonal mean tracer abundance for a particle radius of 1 micron. The 
black lines correspond to the 50% abundance contours. Tracers reach the highest altitudes at mid¬ 
latitude (at around 40°N/S). Except at low pressure for high metallicity (i.e. lOOxsolar and the 
pure water atmosphere), there is a minimum of particle abundance at the equator. As explained 
before, the particle distribution is mostly controlled by the zonal mean circulation with in particular 
the anti-Hadley circulation (see figure]^. On average, particles are thus lofted at mid-latitudes and 
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sediment at the equator. This behavior on a synchronously rotating planet is opposite to that for a 
non-synchronously rotating planet like the Earth, for which there is an upward transport of water 
vapor and chemical species at the ITCZ (Inter-Tropical Convergence Zone), and thus a maximum 
of clouds at the equator. 


Figure 12 shows maps of particle abundance at 12, 1.1 and 0.1 mbar. As discussed in the 
section 4.1, there is almost no longitudinal variation at pressures higher than 1 mbar. The zonal 
wind spatially homogenizes the particle distribution. At pressures lower than 1 mbar, strong 
longitudinal variations appear. These variations are controlled by local upward/downward particle 
fluxes and also by meridional winds, producing more complex distributions. In particular, in the 
upper atmosphere, the meridional circulation tends to transport particles poleward. This is clear 
in figure 12 for H 2 O and the lOOxsolar metallicity case at 0.1-1 mbar. The poleward transport 
is mostly due to the stationary waves, similar to the heat transport in figure Indeed, at low 
pressures, the dayside-nightside contrast is strong and stationnary waves dominate the transport. 
The particle abundance is higher in the dayside, where there is upwelling, than in the nightside, 
where there is downwelling. The poleward wind in the dayside and the equatorward wind in the 
nightside, therefore produce a net poleward particles flux. This meridional transport maintains a 
high particle abundance at high latitudes in the upper atmosphere, and the maximum of particle 
abundance migrates poleward with altitude (see hgure [IT|) . 


4.4. Equivalent ID mixing 

Most studies concerning the photochemistry or the cloud formation on exoplanets are based on 
ID models that take into account the vertical mixing by using a single eddy diffusion coefficient K^z 
to represent the global 3D mixing. In these studies, the values used for the eddy diffusion coefficient 
are quite arbitrary, exploring a range of values and generally without altitude dependence. However, 
our 3D model results can be used to provide eddy diffusion coefficients for use by ID models. 

At steady state, the downward particle flux by sedimentation is compensated by the upward 
vertical mixing. Thus the particle abundance n (in kg/kg or mol/mol) in ID is given by: 


dn 

- PnVf = 0 


( 12 ) 


With the outputs of the GCM, an equivalent ID K^z for the global 3D mixing can therefore 


be defined as (Parmentier et al. 2013): 


{pnVf) 


(pg) 


( 13 ) 


where the brackets correspond to a horizontal and time average over the entire planet. 
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Defining the diffusive time scale and the sedimentation time scale Tg = ^ (both 


depending on altitude), we can also write equation (11) in pressure coordinate as (Parmentier et al. 


2013): 


dn 

'dP 


Tdn 

T,P 


(14) 


If we assume 1) a simple form for the dependence of Kzz on pressure, Kzz=Kzzo{Po/P)°^ with 
a 7 ^ 0 and a (in our simulation a ~ 0.4), 2) an isothermal atmosphere and 3) /3=1.656ifn for 
the terminal velocity in equation (corresponding to a regime where the sedimentation speed is 
proportional to the pressure), then the particle abundance n at a pressure P can be expressed in 


function of no at a pressure Pq as (see also Parmentier et al. (2013)): 


n = noexp 


Tdn 


Tso a-1 



a-l > 
^ 1 -1 


(15) 


This relation can match quite well the mean tracer abundance in figure [TOl for instance for the 
lOOxsolar metallicity using Kzzo = 3 x 10^ m^/s, Pq=1 bar, a=0.4 and no=l. 


Figure 13 shows Kzz obtained from the 3D simulations using equation (12). The left panel 
corresponds to the lOOxsolar metallicity case. Kzz was computed for particle radii from 0.1 to 
3 microns (red lines) and is similar for these different radii. For the dependence on pressure, a 
good parametrization is KzzO = 3 x 10^ x p~^-^ vc?/s. We also compared Kzz obtained from the 
GCM to previous estimates using 1) the simple formula Kzz = WrmsL where Wrms is the root mean 
square of the vertical wind and L is the mixing length assumed to be equal to the scale height H 
( Moses et al.|201~ Lewis et al. poTol ), and 2) the following formula for thermal vertical mixing from 


Gierasch &: Conrath 

(1985 

) and used in the ID cloud model of 

Ackerman &; Marley 

(2001 


Kzz = ^\^. 


H f ( RFc 


\h) \HPaCp) 


(16) 


where R is the perfect gas constant, p is the mean molecular weight, pa is the atmospheric 
density, Cp is the specific heat and Fc is the convective heat flux. The latter is computed from our 
GCM by subtracting the mean net thermal flux by the mean net stellar flux (figure]^. The factor 
1/3 on the left is empirical, derived from observations of giant planets in our own Solar System. 


This formula was also used by Morley et al. (2013) for clouds on GJ1214b 


The equivalent ID Kzz from our model is generally one order of magnitude lower than these 


two estimations. The mixing length has to be lower than the scale height, as predicted by Smith 


(1998). A much better estimate is obtained for both formulae using a mixing length 10 times lower 


than the scale height. Since ID simulations are generally at radiative equilibrium (Tc=0), Morley 


et al. 


(2013) assumed Fc = for the convective heat flux. While such an estimation can be 
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quite good for the convective region of exoplanets or brown dwarfs, it strongly overestimates the 
convective heat flux for irradiated planets. According to our GCM, Fc varies between 1 to 100 
W/m^. The use of Fc = crT^ff overestimates Kzz from formula (15) by one order of magnitude. 


Therefore, the vertical mixing is likely overestimated by one order of magnitude in Morley et al. 


(2013). This will tend to overestimate the particle abundance and size at high altitude in their 
model. Consequently lofting larger particles to higher altitudes to explain GJ1214b’s flat spectrum 
may be more challenging than their results suggest. 


Such simple parametrizations of K^z should therefore be used with caution since they give a 
very rough estimate of the mixing. However, they can give a good trend of the evolution of K^z 
with pressure. In particular, the variations with pressure of Kzz computed with our GCM follow 
quite well those from formula (15). Even if the atmosphere is never convective in our GCM, we 
can expect that the global circulation, driven by the horizontal temperature gradient, globally acts 
as convection for the vertical mixing. The equivalent convective flux Fc in our simulations slowly 
decreases from the altitude where most of the stellar flux is absorbed. For the 100xsolar case, it 

from 0.1 bar to 0.1 mbar. Since pa oc P, equation (15) gives Kzz oc 
) found Kzz oc P~^^'^, while in our model the best fit is with -0.4 as the 
exponent. In fact, an exponent of -0.4 or -1/3 can also provide a good fit to their values. 

The vertical mixing in planetary atmospheres is mostly done by molecular diffusion above 
the homopause and by waves and convection below. On Earth, the vertical mixing is primarily 
done by convection (dry and moist) in the troposphere and by gravity waves in the stratosphere. 
The amplitude of the latter grows as until they break in the mesosphere (Lindzen 1981). 

In warm exoplanets, the mixing should mostly be done by both largescale upwelling/downwelling 
(acting as convection) and planetary waves, which are well simulated by GGMs and generally grow 
as We can therefore expect for such atmospheres an eddy diffusion coefficient with an 

exponent between -1/3 and -1/2 for the pressure dependence, justifying the intermediate exponent 
used in our fit. 


(Lindzen 1981 


is more or less constant 


Parmentier et al. 


(2013 


The right panel in figure [T^ shows the equivalent ID Kzz for the different atmospheric compo¬ 
sitions. It decreases with metallicity in the deeper atmosphere (i.e. at around 1 bar) and increases 
with metallicity above 0.1 bar. This is explained by the difference in altitude where most of the 
stellar energy is deposited (i.e. below 100 mbar for the solar metallicity case and at around 10 
mbar for the 100xsolar metallicity case). The equivalent ID eddy diffusion coefficients can be 
parametrized as Kzz=Kzz 0 x P~^-^ (P in bar) with Kzz0=7 x 10^, 2.8 x 10^, 3 x 10^, 3 x 10^ m^/s 
for the 1, 10, lOOxsolar metallicity and pure water case respectively. The parameters in our fits 
were chosen to primarily match Kzz and the mean tracer abundance (Fig. 10) for pressures lower 
than 10 mbar. At higher pressures, there are quite strong deviations from this simple power-law 
dependence, in particular for the pure water case. For the H-dominated atmosphere, we find that 
Kzz increases with metallicity. A higher metallicity leads to stronger dayside/nightside tempera¬ 
ture contrasts and stronger vertical motions producing a stronger vertical mixing. Kzz for the pure 
water case is lower than for H-rich atmospheres except between 0.1-0.01 bar. However, the pure 
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water case always leads to a higher particle abundance because of its smaller scale height and the 
weaker sedimentation velocity. 


5. Discussion 

The atmospheric circulation and the vertical mixing we described in the last section have strong 
implications for the formation of high clouds/haze on GJ1214b and other tidally locked exoplanets. 
The circulation patterns we obtained, and in particular the anti-Hadley cell, are very general and 
can be applied to a wide range of warm synchronously rotating exoplanets, from mini-Neptunes to 
Jupiter-mass exoplanets. 


5.1. Implications for cloud formation 


The temperature profile in the atmosphere of GJ1214b allows the condensation of potassium 


chloride (KGl) and zinc sulfide (ZnS) (Miller-Ricci Kempton et al. 2012; Morley et al. 2013). Ac¬ 


cording to our simulations for a H-rich atmosphere, the condensation should occur between 40 
mbar and 0.4 bar (see section 3.1.). If the atmospheric circulation can transport cloud particles to 
the upper atmosphere (i.e. 0.1 mbar and less), they should remain condensed (in solid phase) and 


could produce the observed flat transit spectrum (Kreidberg et al. 2014b). Micrometric-size parti¬ 


cles (or bigger) are required to produce a flat spectrum between 1.1 and 1.7 micron. Transporting 
such particles from around 0.1 bar to 0.01 mbar requires a strong circulation that can be accu¬ 
rately quantified only with GCMs. Our 3D simulations with simple tracers can represent quite well 
clouds that form deeper in the atmosphere and do not evaporate at higher altitudes on the dayside. 
They reveal that micrometer particles can be lofted to the upper atmosphere quite efficiently. For 
instance, the abundance is reduced by around a factor 5 at 0.1 mbar for the 100 xsolar metallicity 
with a particle radius of 1 /rm. The vertical mixing tends to increase with metallicity. Therefore, 
metallicity impacts the formation of clouds not only by fixing the amount of condensable species 
but also by controlling the vertical mixing. A high metallicity should therefore favor the formation 


of high clouds on GJ1214b. Assuming the solar nebula atomic abundances from Fodders (2003), 


the abundance of KCl vapor in the deep atmosphere would be around 8.3x10 ®kg/kg for the solar 
metallicity and around 4.3xl0“^kg/kg for the 100 xsolar metallicity. Using the mean profiles of 


tracer relative abundances for particle radii of 1 ^m (Fig. 10), we find that KCl clouds would be 


optically thick in transit at visible and near-infrared wavelength at pressures higher than 3 mbar for 
the solar metallicity and at pressures higher than 0.1 mbar for the 100xsolar metallicity. Therefore 
the latter case is more likely to produce the observed flat transit spectrum than the former case. 


Clouds could also modify the circulation and the vertical mixing by absorbing and scattering 
stellar radiation. If clouds are strongly scattering, we can expect a high planetary albedo and 
in particular a cooling of the atmosphere below the cloud deck (i.e. below 0.1-1 bar) impacting 
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the atmospheric circulation. The impact of this cloud-climate feedback on the possibility of high 
clouds on GJ1214b is best understood using a GCM that can take into account the radiative effect 
of clouds. Finally, our GGM predicts a minimum of tracers at the equator for pressures higher 
than 0.1 mbar for all cases considered here, and we expect the cloud distribution to have the same 
behavior. A minimum of cloud at the equator could have observable impacts on transmission 
spectra and also on reflection or emission spectra of the planet. 


5.2. Implications for photochemical haze formation 


Another possible explanation for the flat transit spectrum of GJ1214b is the formation of 
photochemical haze, possibly produced by Titan-like methane photolysis. Such a process can be 
more efficient than cloud condensation at producing high aerosols because particles are produced 
directly in the upper atmosphere. However, sedimentation may be too strong to maintain a particle 
density large enough for blocking stellar light. Even on Titan, which is very hazy and has a 
low gravity (with correspondingly small sedimentation velocities), the haze is not opaque enough 


to block all methane absorption bands in near-infrared transit spectra (Robinson et al. 2014). 


Therefore GJ1214b’s flat spectrum can only be produced by photochemical haze if there is a large 
methane flux counter-balancing the haze sedimentation flux. Based on the atmospheric circulation 
simulations, we predict that methane would be supplied at mid-latitudes (see figure [^. Once 
produced, haze would be removed mostly at the poles and at the equator. This would impact the 
hazes’ latitudinal distribution. Haze particles should fall until the temperature is high enough to 
pyrolyze them, mostly into methane. The Huygens Probe Aerosol Gollector and Pyrolyser (AGP) 
experiment showed that Titan’s hydrocarbon haze was completely pyrolized at temperatures close to 
900K ( Israel et al.||2005 ). In our simulations of GJ1214b, that temperature corresponds to pressure 
of around 0.1 bar for high metallicity and 1 bar for solar metallicity (see Fig. [^. Therefore, there 
could be a complete methane cycle for pressures lower than these values, that could recycle haze into 


methane. Miller-Ricci Kempton et al. (2012) and Morley et al. (2013) explored the photochemistry 
in GJ1214b with a ID model using eddy diffusion coefficient values from 10® to 10® cm^/s and 
predicted that more methane and more haze would be present at high altitude with a high K^z 
(i.e. 10® cm^/s). With our model, we predict eddy diffusion coefficient values of around 3x10^ 
cm^/s at 1 bar to 8x10® cm^/s at 10“® bar for the lOOxsolar metallicity case. The photolysis of 
methane, which occurs mostly above 10“® bar, is thus better simulated with A' 2 ^= 10 ® cm^/s in 


Miller-Ricci Kempton et al. (2012) and Morley et al. (2013). Moreover the radiative effect of haze. 


which can absorb stellar radiation, could increase the vertical mixing in the upper atmosphere and 
support the lofting of aerosol particles, increasing the lifetime of the haze, and this feedback is best 
studied with a GGM. We conclude that the formation of high haze on GJ1214b is likely and it is 
also possible that both high condensate clouds and photochemical haze are formed. 
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6. Conclusions 


In this study, we analyzed the atmospheric dynamics of GJ1214b for a cloud-free atmosphere 
with different metallicities as a first application of the Generic LMDZ GGM to a non-rocky exo¬ 
planet. We obtained results for radiative transfer, temperatures and winds very similar to those 
from the other state-of-the-art GGM (Kataria et al. 2014), validating our model for this kind of 
planet. 


In addition to this intercomparison, we showed that the zonal mean circulation of GJ1214 
is characterized by the presence of an anti-Hadley circulation on the nightside, which leads to 
zonally-averaged subsidence at the equator. This particular regime likely occurs on a wide range 
of synchronously rotating exoplanets, from warm mini-Neptunes to hot Jupiters. 


Using simple tracers, we analyzed the vertical mixing in the atmosphere of GJ1214b. The 
tracers showed that atmospheric transport is primarily driven by the zonal mean circulation, leading 
to an upwelling at mid-latitudes and a downwelling at the equator. In particular, a minimum of 
tracer abundance appears at the equator, strengthened in the upper atmosphere by a poleward 
meridional transport. This should be a fundamental feature of the cloud/haze distribution on many 
synchronously rotating warm exoplanets. We also found that the vertical mixing increases with 
metallicity for H-rich atmospheres. For a solar metallicity, the circulation cannot loft micrometer 
particles into the upper atmosphere. Therefore, if the upper atmosphere of GJ1214b is opaque 
because of condensate clouds, the atmospheric metallicity must be higher than solar and likely 
more than 10 xsolar. From these simulations, we derived equivalent ID eddy coefficients. We found 


- 0.4 


that: Kzz = 7 x 10^ x P^^.' m^/s is a good fit for lx solar metallicity and Kzz = 3 x 10^ x 
m^/s for 100 Xsolar metallicity. We compared these values to classical formulae for Kzz used for 
instance in Morley et al. (2013), showing how these simple formulae can overestimate the mixing 
by one order of magnitude. Our parametrizations of Kzz can be used in ID cloud or photochemical 
models ( Miller-Ricci Kempton et al.||201^ Morley et al.|[2M3 ) for which they are expected to lead 
to an efficient formation of photochemical haze. 


Finally, we showed that the atmospheric circulation of GJ1214b could be strong enough to loft 
micrometric particles of KGl or ZnS from 1 bar to 0.1-0.01 mbar. However, the radiative effects of 
clouds, by absorbing or diffusing stellar radiation, could modify the strength of the circulation. The 
next step will be to simulate the atmosphere of GJ1214b with realistic KCl and ZnS clouds, taking 
into account the latent heat release and the radiative effects. For this next study, the Generic LMDZ 
GGM benefits from all the developments done for previous studies of cloudy rocky planets. Such 
3D simulations would give strong indications of the possibility of high condensate clouds producing 
the observed flat transit spectrum. Such simulations could also help identify the best observations 
for probing cloudy atmospheres, with, for example, emission/reflection spectra and phase curves. 


B.G. acknowledges support from an appointment to the NASA Postdoctoral Program at NAI 
Virtual Planetary Laboratory, administered by Oak Ridge Affiliated Universities. VSM is sup- 











19 - 


ported by the NASA Astrobiology Institute’s Virtual Planetary Laboratory Lead Team, under 
the National Aeronautics and Space Administration solicitation NNH12ZDA002C and Cooperative 
Agreement Number NNA13AA93A. This work was facilitated though the use of advanced compu¬ 
tational, storage, and networking infrastructure provided by the Hyak supercomputer system at 
the University of Washington. 



- 20 - 


Parameters 

Dp (planetary rotation rate, rad s“^) 

4.615x10“^ 

Prev (revolution period, s) 

136512 

e (eccentricity) 

0 

Rp (planetary radius, m) 

1.7x10^ 

g (gravitational acceleration, m s“^) 

8.93 

Fstar (stellar flux at top of the atmosphere, W m“^) 

23600 

Pint (internal thermal flux, W m“^) 

0.73 


Table 1: Orbital and physical parameters used in the model. 


Atmospheric composition 

Cp (J kg-i K-^) 

L7(km) 

Mean molecular weight (g mol ^) 

1X solar 

12800 

220 

2.3 

10 X solar 

12000 

200 

2.5 

100 X solar 

6474 

115 

4.38 

pure water 

1864 

28 

18.0 


Table 2: Values of specific heat (cp), scale height (H) and mean molecular weight used for the 
different atmopheric compositions. 



Fig. 1.— Left panel: temperature profiles from the ID model for the different atmospheric compo¬ 
sitions. Dashed and dotted lines correspond to the saturation vapor pressure curves for KCl and 
ZnS, for the solar (gray) and the 100xsolar (black) metallicity. Right panel: maximum dayside- 
nightside temperature difference as a function of pressure. The differences are computed at each 
latitude and pressure level. The figure shows their maximum value at a given pressure. 
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Fig. 2.— Global mean heating/cooling in K/s (left) and stellar/thermal net flux in W/m^ (right) 
as a function of pressure. Solid lines correspond to the heating rate and the stellar net flux. Dashed 
lines correspond to the cooling rate and the thermal net flux. 
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Fig. 3.— Zonally averaged zonal wind in m/s for the four atmospheric compositions. 
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Equator-Temperature (lOOxsolar) Equator-Vertical velocity (lOOxsolar) 



Longitude (°) Longitude (°) 


Fig. 4.— Temperature (in K) and vertical velocity (in m/s, positive value = upward wind) for the 
lOOxsolar metallicity. The top panels show the temperature (left) and vertical velocity (right) at 
the equator versus longitude and pressure. The bottom panels show the map of temperature (left) 
and vertical velocity (right) at 1.1 mbar. Black vectors correspond to the directions of horizontal 
winds. 
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Latitude (°) 

Mass streamfunction (nightside, tOOxsolar) 



Latitude (°) 


Fig. 5.— Zonal mean mass streamfunctions in kg/m/s for the lOOxsolar metallicity case for the 
dayside (top), the nightside (middle) and globally (bottom). Positive (negative) values correspond 
to anti-clockwise (clockwise) circulation. Dashed (solid) lines are contours for the value of ±10^^ 
kg/m/s. The dayside streamfunction was calculated by integrating longitudes between -39° and 
-|-90°E. The nightside streamfunction was calculated by integrating the other longitudes (-180 ° 
to -39°E and -|-90° to -|-180° E). The global mass streamfunction was calculated by integrating all 
longitudes and is equal to the sum of the dayside and nightside streamfunctions. 
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Fig. 6.— Zonally-averaged vertical wind in m/s (positive value = upward wind). Solid (dashed) 
lines correspond to contours of the zonally-averaged mean mass streamfunction with clockwise 
(anti-clockwise) rotation. 
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Fig. 7.— Heat transport by the meridional atmospheric circulation. Top panel: zonal mean differ¬ 
ence between the outgoing emitted flux and the absorbed stellar flux for the different atmospheric 
compositions. Positive (negative) values correspond to regions globally warmed (cooled) by the 
circulation. Bottom panel: decomposition of the total heat transport (black solid line) for the 
100 X solar metallicity case into heat transport by the mean meridional circulation (black dashed 
line), by stationary waves (black dashed-dotted line) and by transient perturbations (black dotted 
line). Positive (negative) values correspond to northward (southward) heat fluxes. 














- 26 ~ 




Fig. 8.— Variations with pressure of the sedimentation velocity (top) and the sedimentation rate 
(bottom) for particles with radii from 0.1 to 10 microns. The solid and dashed lines use the ID 
temperature profiles for the 100xsolar metallicity (solid) and the pure water atmosphere (dashed). 
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Fig. 9.^— Zonally-averaged ratio of the sedimentation timescale by the advection timescale for the 
100 X solar metallicity. Contour colors are in log scale. 
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Tracer abundance (r=1 ^im) 



Tracer abundance (100x solar) 



Fig. 10.— Mean profiles of the relative abundance of tracers. Left panel: for the different atmo¬ 
spheric compositions with particle radius of 1 micron. Right panel: for the 100 x solar metallicity 
with particle radii from 0.1 to 10 microns. For both figures, the relative abundance is fixed to 1 
below 1 bar. The dashed lines correspond to the analytical profiles using a parametrized Kzz (see 
section 4.4). 
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Fig. 11.— Zonally-averaged relative abundance of tracers with particle radius of 1 micron for the 
different atmospheric compositions. The black line corresponds to pressures where the relative 
abundance is 50%. 
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Longitude (°) 

lOOxsolar, 0.1 mbar, r=1 [im 



Fig. 12.— Maps of relative tracer abundance (particle radius of 1 micron) for the different atmo¬ 
spheric compositions (lx, lOx, lOOxsolar and pure H 2 O from top to down) and for pressure of 12 
mbar (left) 1.1 mbar (middle) and 0.1 mbar (right). 
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Fig. 13.— Profiles of equivalent ID eddy diffusion coefficients. The left panel shows for the 
100 Xsolar metallicity case. Red lines are from the GCM for particle radii of 0.1, 0.3, 1, 3 microns. 
The bine line is the simple fit {K^z = 3 x 10^ m^/s). Solid and dashed black lines 

are estimations from simple formula (jGierasch &: Conrath 1985[ Ackerman &; Mar ley 2001; Lewis 


et al. [2010 Moses et al.||2011 ). The right panel shows ¥^zz derived from the GCM for the different 
atmospheric compositions with particle radins of 1 micron. The dashed lines are the simple fits 
{Kzz=Kzz 0 X with Kzzo=7 x 10^, 2.8 x 10^, 3 x 10^, 3 x 10^ m^/s for the 1, 10, lOOxsolar 

metallicity and pure water case respectively). 
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